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Abstract 

In this paper we construct a pseudorandom multisequence (x nii ,„ iTh .) 
based on A;th-order linear recurrences modulo p, such that the discrep- 
ancy of the s-dimensional multisequence (xni+u,...,n r +i r )i<i J <s J ,i<j<r 
1 < nj < Nj, 1 < j < r is equal to 0((Aq • • • N r )- 1 / 2 ln s+3r '(iVi • • • N r )), 
where s = s\ • • • s r , for all N\, ...,N r with 1 < Ni ■ ■ ■ N r < p k . 



1 Introduction 



Equidistribution and statistical independence properties of uniform pseudorandom 
numbers can be analyzed based on the discrepancy of certain point sets in [0, l) s : 

Let x n = (x nt i, . . . , x ntS ), n = 0, . . . , N — 1, be a sequence of points in an s- 
dimensional unit cube [0, v = [0, 71) x . . . x [0, 7 S ) a box in [0, l) s . The quantity: 



sup 

0<7i,...,7 s <l 



#{ne [0,JV-1] I x n e^}/iV-7 1 ...x t 



,N-1 



is called the discrepancy of (x n ) n=0 . 

Let us consider pseudorandom numbers (abreviated PRN) obtained by means of 
the shift-register method: 

Let p be a prime, let k > 2 be an integer, and generate a fcth-order linear 
recurring sequence yo, yi, ... G F p by 



y n+k = a k _iy n+k _i + ... + a y n mod p, 



n = 0,1,..., 



(1) 
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where yo>--->yfc-i are initial values that are not all zero. The integer coefficients 
do, ...,Ofc_i in ([1]) are chosen in such a way that, if they are viewed as elements of 
the finite field F p , then the characteristic polynomial 

f(x) = x k — a k _i% k ~ l — ... — a E F p [x] 

of the recursion (JTJ is a primitive polynomial over F p . Note that the characteristic 
polynomial / has a root /3 in the extension field F q of F p , where q = p k . Let 
F* — F q \ {0} be the multiplicative group of nonzero elements of F g and let F q = 
{f3 e F* | — is a primitive root}. We see that #F* = q — 1, and #F 9 = 
(^(g — 1), where ip is the Euler's function. Let Tr denote the trace function from F q 
to F p . It is known (see, e.g., [Ni, p. 212]) that there exists an a G F* such that 

y n = Tr{al3 n ) for n = 0, 1, ... (2) 

In the digital multistep method, the sequence yo,yi, ■■■ is transformed into a se- 
quence Xo,Xi, ... of uniform PRN in the following way 

k 

X n = ^Vkn+j-l/p 1 , for 71 = 0,1,... (3) 

In a series of papers, Niederreiter (see the review in [Ni]) proved that there exists 
a characteristic polynomial / such that 

D((x n ,...,x n+s . 1 ^- 1 ) = O( ^ iln J r+1 ), for N = l,...,r, 
where r is the period of the sequence of pseudorandom numbers. 

This estimate is interesting for iV ^> -^(log r) s+1 . In [Lel],[Le2], Levin descibed 
a class of uniform PRN sequences (z n ) n > , having a nontrivial discrepancy estimates 
also for a small part of the period: 

D((z n ,...,z n+s . 1 )^- 1 )=O(N- 1 / 2 ln s+3 N), for iV = l,2,... (4) 
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Our goal is to obtain a nontrivial discrepancy estimate similar to (jl]) for a small 
part of the period for sequences, some subsequences and multisequences of PRN 
based on fcth-order linear recurrences modulo p. The method of the proof is based 
on Korobov's approach [Ko2] (see also [Lei], [Le2]). Similar results can be obtained 
for the pseudorandom sequences described in [Le2] and [Le3]. 

In this paper we will prove the following theorems: 

Theorem 1. Let e G (0, 1), and let x(n) = (x n , . . . , x n+s _i) (see ^j, ^)). Then 
there exist more than (1 — e)q(p(q — 1) pairs (a, (3) G F q x F q such that for any 
N G [l,q] and any 61,62 £ [1>9]> ^e bound 

D((x(6 1 n + 6 2 ))£r 1 ) < e^dN-^Hn^ 2 - 5 6N\n 2 - 5 3\n6N, 

holds, where the constant c\ depends only on s, 6 X and 6 2 . 



Let r > 2, a = (ai, . . . , a T ) G F£, and (3 = ((3 1 , . . . , /3 r ) G FJ" , n = (ni, . . . , n r ) 
and let 

r 

^j^^Tr^"^^- 1 ) modp, y(n,j)e[0,p-l], 

i=l 

where, 

n.j = n w s dwa mod s r_1 , n w G [0, s r_1 ), 

iw£[l,r],w^i 

with 

d w ,i — w — 1 for w < i, and = w — 2 for w > i. (5) 
Consider the following multisequence 

— , and x n = (x n+i )o<i„< Slu , with S; G [1, s], 1 < w < r, (6) 

. , p 3 Kw<r 
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where i = (ii, . . . , i r ). 



Theorem 2. Let e G (0, 1). Then there exist more than (1 — e)q r ((p(q — l)) r 
values (a, (3) G F r q x F r q such that, for any N { e [l,q], 1 < i < r, N x . . . N r < q, 
and any s$ G [1, s], 1 < i < r, s = si . . . s r < s, the bound 

D((x n )o< nw <N w ) < e^c 2 {N ± . . . N r )~ l/2 \n So+2 - 5r (2 r+1 N, . . . N r ) ln 2 ' 5r 3 ln^ • • • N r ) 

l<w<r 

holds, where the constant c 2 depends only on s±, s r . 

2 Auxiliary results 

For the integer b > 2, let denote C(b) = (-b/2, b/2] f| Z. Let C s (b) be the inner 
product of s copies of C(b). Consider point sets for which all coordinates of all 
points have a finite digit expansion in a fixed base b > 2. Let 

w n = («,«, . . . , G [0, l) s , n = 0, 1, . . . , N - 1, (7) 

where, for an integer m > 1, we have 

m 

wf = ^ w^b~ j , < n < N - 1, 1 < i < s, 



define 



with w% G Z b for < ra < N-l, I < i < s, 1 < j < m. For . . . , h m ) G C m (6), 



d(hi,...,h r 



largest d with hj^O, if (/ii, . . . , h m ) ^ 0, 
0, if (/n,...,/i m ) = 0. 



For 6 = 2, put ...,h m ) = 2-< hl '-' hm \ For 6 > 2, put 

Qb(hi, . . . , /i m ) = 



6- d (cscf|h d | +<7(d,m)), if (hi,...,h m )^0, 
1, if (h 1: . . . , h m ) = 0, 
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where d = d(h\, . . . , h m ) and where a(d, m) = 1 for d < m and a(m, m) = 0. 

Let C(b) sxm be the set of all s x m matrices with entries in C(b), C*(b) sxm = 

C(b) sxm \ {0}. For H = {h i3 ) e C(b) sxm , we define 

s 

W b (H) = Y[Q b (h il ,...,h im ). 



Theorem A. [Ni, Theorem 3.12, p. 37, Lemma 4.32, p. 68] If P is the point set 
Ffy, and m > [log 6 N] then 



sb 



N-l ... 

^E e (iEE ft 



n=0 



i=l j=l 



where e(x) 



^2nV-lx 



Lemma 1. [Ni, Lemma 5, p. 18] Let s > 1 and m > 1 be integers. Then 

Y W b (H) < (-m\nb+-m-^^) . (8) 

^— ' V 7T 5 b 



HeC*(b)° 



Lemma 2. (see e.g., [KoSh, p. 9, p. 13, ref. 3.3]) Let f3 G F q , 

6{J3) 

Then 



1, if /3 = 0, 
0, otherwise. 



Tr(«/3) 
P 



For proof of the following well known lemma see, e.g., [Ko, p. 13], or [LeVo, Lemma 
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7, p.156]. 

Lemma 3. Let N G [0, T — 1], T G [1, g] and x ra fre a rea/ (t) < n < T — 1). Then 

T/2 



7V-1 



ra=0 



< V - Ve(i„ + 
^— ' m z — ' V 

m=-T/2 n=l 



nm\ 
~f~) 



where m = max(l, \m\). It is easy to see that 

T/2 



V — <3 + 21nT. 



m=-T/2 



(9) 



Lemma 4. Lei G [0,1* - 1], 1 < i < r ; w/iere T; G [l,g], 
andn= (ni,...,n r ). JTien 



iVi-l 7V r -l 

X) ■ ■ ■ Z) 

ni=0 n r =0 



Ti/2 

S E 



T-/2 

E 



X 



mi=-Ti/2 m r =-T r /2 
Z • • • Z 6 ( Xn + 



mi . . . m r 



Ti-l T--1 



nimi n r m r 
+ " ' + T r 



T-i 



ni=0 n r =0 

The proof of the Lemma 4 is the same as the proof of the Lemma 3. 



1 < i < r, 



Lemma 5. Let r > 2, S\, . . . , s r G [1, s] be integers, 

m+ii= ^2 (n w + i w )s dw - 1 mod s r ~\ n t + i x G [0, s r_1 ) 



wS[l,r] 



Then for I G [1, r], 



(10) 



< j < < < s^, v = I, . . . ,r j = ks\ . . . s r . 

Proof. By ffTOl) it is enough to prove that there are no two vectors (ii, . . . ,i r ,j) ^ 
(ii,...,i' r ,j') with 



k(sni + ii) + j = k(sni + i{ + i[) + j' . 



(11) 



We see j = j mod k. Hence j = j' . By (TTTj) . we get: 

i\ = i t mod s. (12) 
Therefore %\ = i[. From ffTTj) we have that n t + %i — ni + i\. Hence 



( n ™ + ^) S<W = Y (n w + i' w )s dw - 1 mod s 



r-l 



Thus 

{i w -i'js^' mod^ 1 . (13) 



tuS[l,r] 



For I 7^ 1, we obtain 

— z : = mod s. (14) 

Note that for I — 1, (TJJJ follows from ffl2l . Hence = i x . Suppose that ij = i'j for 
j = 1, . . . , v — 1, v ^ Z, i/ > 2. By ([5]) we get g^z — o^z > 1 for u> > z/. We deduce 
from ([TJ) that 



v - C + Y ~ i '^ sdwJ du ' 1 = mod s 



r—l—du 



wE(v,r ] 



i v — i u = mod s„, and % v = i v . 



By induction, Lemma 5 is proved. □ 

Lemma 6. Let q = p k > 3000, T G [l,4g]. Then 

< 401n3Tln(31n3T). (15) 



<p(q-l) 
Proof. First we will prove that 

T 



. < 10 In 3 In 3T. (16) 
1) 
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By [Sa, p. 15, ref. 3a; p. 9, ref. 2] 



n 2 6 

< e - 58 In Inn + — — , and <p(n) > n 2/3 n > 30. 



ip(ri) In Inn 

Hence 



<-j + < 1.5 + 2 lnln(g- 1). (17) 



<p{q-l) ~ <p(q-l) (g-1) 2 / 3 
If T < <p(q — 1), then (fT6l) is true. Let T > ip(q — 1), then 

2 

InT > -ln(g- 1), and lnlnT > lnln(q - 1) - 0.5. 

3 

By (fT7|) we have 

T 



<4(1.5 + 21nln(g-l)) < 10 + 8 In In 3T. 



if(q-l) 

The inequality (fl~6l) is proved. If T < y?(g — l)/fc, then (fl5]) is true. Let T > 
-l)/k> <p(q - l)/log 2 g, then 

InT > lng — ln(1.5 + 21nln(g — 1)) — ln(log 2 g). 

Hence 

41nT > 9(q), where 0(x) = 41nx - 41n(1.5 + 21nln(x - 1)) - 41n(log 2 x). 

It is easy to verify that (9(x) - log 2 (x))' < for x > 3000, and that 0(3000) > 15 > 
log 2 3000. Thus 

4 InT > log 2 g. 

Applying ffl6|) we get ( Tl5|) . Lemma 6 is proved. □ 

3 Proof of Theorem 1 

By (jSJ) and Theorem A, with m = [logpiV] we get: 

N-D((x(b 1 n + b 2 ))^ 1 )<sp+ W P (E)\S(H)\, (18) 

HeC*(b) sxm 
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where 



N-l /_, a— 1 m 



Using Lemma 3, we have 

*(//) < £ — S(H,v ») 



n=0 V i=0 3=1 



T/2 



mi=-T/2 



mi 



where 



with 



S(if, m 1; a, 0) = J2 e (-Tr (a/3 fc ^ 7 (/3, H)) + ^f-) , 

n=0 ^ ' 



s— 1 m 

70S, #) = E E M" 1 *" 8 "" -1 . and T e [N, q]. 

i=0 3=1 

Taking m = [log„T], we obtain from ([15]) and (JT9]1 . that 



N ■ D((x(6m + 62)^=0) < + TZ} T («, /3), 



where, 



T/2 

TD T (a,f3)= W p( H ) E — S(H,m,a,P) 



,ffec*(ft) s 



mi=-T/2 



Let, 



.T^pryEEl^™'^) 



Using the Cauchy - Shwartz inequality, we get: 

T-l 

V V 

q ■ <p(q 



x 2 < 



, _ 1) E E E 



T-l 



(19) 



(20) 



(21) 



(22) 



Tr (a/3 fcbini 7 (/3, #)) - Tr (a/3 fcbin2 7 (/3, //)) ' 



/3G-F, a ^ F g ni,n 2 =0 



P 



xe 



mi(ni - n 2 , 
f 



Applying Lemma 2, we obtain: 

T-l 



X 



<p(q 



3T)E E 5(7(/3,^)(/3 fc6ini -/3 fcfein2 ))e( 



iSg ^ni,n 2 =0 



mi(ni - n 2 , 
f 



T-l 



Bearing in mind that 5(7172) < S(ji) + £(72), we have: 



X < Xi + X2, 



where 



Xi 



<p(q 



1 T_1 

31) E E t(P kbini -P kbin2 ), and X 2 



/3eF 9 ni,n 2 =0 



with 



(23) 



(24) 



Consider xi- Let f3 kbini - f3 kbin2 = 0. Then /3 fc M«i-™2) = L Taking into account 
that (3 ia a primitive root, we obtain kb\{n\ — n 2 ) = mod y>(g — 1). Hence, 



Xi < T 1 + 



(25) 



Now consider % 2 - By ([HI . (|2"3|) . £ is equal to the number of solution of the following 
equation: 

s— 1 m 

(26) 



i=0 j=l 

Bearing in mind that (/in, . . . , /i sm ) 7^ 0, and (126]) is a polynomial equation on the 
field, we get: £ < fcfcis. Thus 

kb lS T 2 



X2 



< 



<p(q-lY 



(27) 
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By (J23D, (J2SD, (J2ZD and Lemma 6, we have 

kb x T 



r < t i + 



kb lS T 2 m / , , kb x T 
+ / - <T (l + ( a + l). 



< T(l + (s+ l)6i401n3Tln(31n3T)) < 40(s + 2)&iTln3Tm(3 ln3T). 
We see that 

2 7 m — 1 

— mlnpH — m < 2.5 InT, for m = [log„T] > 1. (28) 

7T 5 p p 

From (I2"T1) . (|22|) . ((HD, ©, Lemma 1 and Lemma3, we obtain: 

— r V V TD T (a,f3) < (2.51nT) s (3 + 21nT)40 1/2 (s + 2) 1 / 2 feJ /2 T 1 / 2 

qip(q — 1) z — ' ^— f 

xln 1/2 3Tln 1/2 (31n3T) < 3-40 1/2 (s+2) 1/2 2.5^ /2 T 1/2 ln s+L5 3Tln 1/2 (3 ln3T). (29) 



Let Tj = 4*, i > [log 4 p] + 1, and let 

q q [log 4 4q] 



J ^ £l i=[1 ^, ]+1 40V2( S + 2)1/22.5- ln s+2 ' 5 3T, ln 25 (3 In 3T,) 



By pgp, we have: 



q-tp(q- 1) 



<J q oo 



< 1. 



Let 

Q e = {aeF q ,peF q 

Let's prove, that 



, H f-T' 96i6 2 ln 2 36iln 2 36 2 3iln4 In 2 (i In 4) 

6 1 = 1 &2=1 i=l v ' 

' i?(a,/3)<l/e}, #Q e = 7 g-^( g -l). 
7 > 1 - e. 



(30) 



(31) 

(32) 
(33) 
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We see that q-<f(q— 1)(1— 7) is the number of a G F q , (3 G F g , such that R(a, (3) > -. 
From (l3Tj) we obtain 



:(i -7)9- 1) = - — -■ 



q-<p(q- 1) e 
The inequality f )33|) is proved. 

Now, let iV G [Tj ,T io+1 ) for some i G [log 4 p, log 4 q], where T t = 4\ q > 3000. 
From (I20D, (EUD and (ETJ we have for all (a,p)eSl e , 

ND{x{n, a, P)) < n <N-i < sp + T l0+1 D TtQ+1 (a, /3)<sp + e'^Ho 1 ^ + 2)^ 2 2.5 S 

x T J 1 / + 2 1 ln s+2 - 5 3T i0+1 ln 2 ' 5 (31n3T l0+1 ) < e _1 ciiV 1/2 ln s+2 ' 5 6iVln 2 - 5 (3 ln6iV), (34) 
with 

c x = sp + 3 6 (s + 2) 1 / 2 2.5 s b 1 1 % In 2 3&i In 2 36 2 . 

It is easy to see that c\ > max(p, 3000). Hence, if A" < max(p, 3000), then ( 134"1) is 
also true. Theorem 1 is proved. □ 

4 Proof of the Theorem 2 

By §6§ and Theorem A, with m = [log p Ni . . . N r ] and s = si • • • s r , we get: 

l<n w <N w , l<w<r )<s oP + W P (H)\S(H)\, 

-ffec*(p) s o xm 

where 

/ , \ 

S(H) = J2 e 

n<e[0,JVi-l] 



k r 

,k(ni+ii+sni+ii)+j-l 
P • ■ 1 



- j . iwe[o,s m -i] j=l 1=1 . 

X,...,r \ «i=l,...,r / 
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Let Ni 6 [l,Tj], Tj <G [l,g], 1 < i < r. Using Lemma 4, we get: 



Ti/2 

w< E 



^V 2 |£(tf,T,«,/3) 



m 1 =-T 1 /2 m r =-T r /2 



mi . . . m r 



where 



with 



S(H,T,a,0) 



ti=0 t r =0 



(35) 



Sl— 1 s r — 1 fc 



^(t, «, /3) = E (- E • • • E E (m a 

!=1 " ii=Q i r =0 j=l 

Hence, 

iVi . . . N r ■ D((x n ) < nw<Nwt i< w < r ) < s p + 7i . . . T r D Tu ... tTr (at, (3), (36) 



where, 



Ti/2 



T 1 ...T r D Tl ,..., Ti .(«,/3)= E E 



^ |£(F,T,<*,/3) 



H&C*(p) a o> 

with m = [log p (Ti . . . T r )]. Let 

X(T) 1 



q r (cp(q-l)y 



mi=-Ti/2 m r =-T r /2 



E E |5(^,T,a,/3) 



mi . . . m r 



(37) 



Using the Cauchy - Shwartz inequality, we get: 

Ti-l T r -1 



X 2 (T) < 



1 



<r • - i)y 



EE 



^...^e(a(t,a,/3)) 



By (|35|) . we have 

X 2 (T) 



t 1= o t r =o 



Ti-l 



(38) 



Tr-l 



<f ■ (<p(q ~ 1))' 



EE E ••• E 1 



,3eFrCXEFr t M jt m =0 t (l) )t (2) =0 
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xe 



\ E ■ ■ ■ £ £ E k > ( Tr ( ai/J <w<> + ,« V 



i 1= i r =0 j=l 1=1 

xe ( -Tr [ ai ^^^n) + j-^ + m^- t^) + + ^M") 
Using Lemma 2, we get 

, Ti-l T r -1 

*>^i: e - e 

jj^x; 1 y~jy^^- ^Hti 2) +n+st\ 2) +ii)+j-i _ pkitw+ii+st^+iti+j-i 

i=i n=o j r =o j=i 

xe I — h . . . + 



I) T r 

Ti-l T r -1 

- 1))' ' 



< • • _w E E ••• E 



/3eF, r t«tf»=o 4 1) ,4 2) =o 
II 5 ( E • • • E E ^ ' (/3* (t « (a) ^ + ^ ,)+i - 1 - ^W^O+.-i^ 

1 = 1 ii=0 ir=0 J=l 

It is easy to see that 

e iW 

«=l.-,ri=l,...,r, j=l,2 



with 

EV^V^A ( p k (ti 2) +st\ 2) +ii)+3-l Mt^+k+st^+i^+j-l 

■ ■ ■ 2_- ''" [Pi ~~ & 

il=0 i r =0 3=1 \ 

We take a new variable v t instead of tf^ (I — 1, ... ,r). Enlarging the domain 

of the summation, we obtain: 

* 2 m s j^zjw E E II E m. 

»=l.-,i-i=l,...,r, j=l,2 
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where 

jkffi+svt+iti+j-l ofc(tj (1) +ii+st t C1) +ii)+i-l 



& = E • • • E E ( # - # (t; ) . (39) 

«1=0 i r =0 j = l 



Hence, 



X 2 (T)<II E 6(t (1) ), (40) 



i=1 tf ^[O.Ti-l] 



with 



«t (I) ) = ^£ E E '<«■ («) 

Ae^tfe^-il^efo.^- 1 ) 



Consider the equation £j = 0. By fl39l) . we have: 



A^-^^i = 72, (42) 



where 

a\— 1 s r — 1 A; 



71 = E • • • E E h h,.,ir,jPl 
il=0 i r =0 j = l 

si— 1 s r — 1 fe 



fe(s«i+ii)+jf-l 



7 2 = 72(A) = E • • • E E ^,..,v.-/r t! (43) 

ii=0 i r =0 j=l 

Similarly to (|23|) . we get: 

6(t (1) )<6(t (1) ) + 6(t«), (44) 

where 

«' (1) ) = ^E E E *(a* <( '' p> ~''" >) 7i-72 

x(l-5( 72 )), and £(t«) = Ti— £ £ 5( 72 ). (45) 

Consider £,(tW). We see that if 72 = 0, then &(t«) = 0. By (|32J) and gSD if 7 2 ^ 
and 71 = 0, then also ^(t^) = 0. Now let 71 7^ and 72 7^ 0. There exists an 
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integer a, such that j3? = 72/71- By fflSl) . we get /c(t| 2 ^ — t[ ) =a mod 
Hence 



#|0<t{ 2) <T ; fc(4 2) -4 1} ) = a mod^(g-l)} < 1 + 
By (Ho^ and Lemma 6, we get 

6 < s r_1 (l + 401n37]ln3(ln37})). 
Similarly to f|27|) . we have 

,2r-l T * 



# {A e F,| 72 (/3,) = 0} < £;s r , and 6(t«) < A; S 2 



p(g-l) 

By f|40l) . (l46l) and Lemma 6 we obtain 

&(t (1) ) < s"" 1 + AOs^ 1 In 37} ln(3 In 3T,) + 40s 2r " 1 In 37} ln(3 In 3T,) 
< 81s 2r - 1 ln3T i ln(31n3T i ). 

From ( l4*0l) . we have 



X 2 (T) < J];T i 81 S 2r - 1 ln3^1 n (31n3rO. 

i=l 

Using (1281) . (1371) . fl38|) . Lemmal and Lemma4, we deduce 

^TTW E J2 T ^-- T r f) T 1 ,...,T r (^/3)<(2.5\n(T 1 ...T r )Y 



q r ((f(q - l)Y 



x JJ (2 + 3 ml)) ^^(Sls 2 ^ 1 hi37] ln(3 ln37})) 1/2 for T 1 ---T r >p. 

1=1 

Let Tj 4 = 4 Ji , ji = 0, 1, i = 1, ...,r, and let 



!< Sl tr Sr < s 1^,..^, 3- Sl ... Sr ln 2 3 Sl ...ln 2 3 Sr ' 

Si — s r <s log 4 p<ji + . . .+Jr <log 4 Aq 
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where 

r 

V>(Tj) = (2.51n(Ti...T r )) s fJs r - 1 / 2 ln 2 - 5 3T ii In 25 (3^3^). 

1=1 

From (T4"T|) we get 



r oo oo 

< " 



nyy — a — ^ — ,<i. 

^ ^ 3s, In 2 3s 4 • 3* In 4 ln 2 ^ In 4) 



Let 



i?(«,/3)<-|>, #a = 7? r M?-i)) r - (19) 



Similarly to (!31~!) - (!33l) . we get 7 > 1 — e. 

Now, let Ni G [T Ji; Tj 4+1 ) for some ji G [0,log 4 g], z = 1, ...,r, with 2^ = 4 Ji and 
Tjj . . . T jr > p. From flU, gHD and (|3SD, we have for all (a, /3) G 

) < s p + . . . T jr+ iD T . (a, /3) 

< SoJ9 + e -l 3 5r a r»-r/2 (^...^1/2 
r 

(2.5 In (T n+1 . . . T jr+l )) s ln 2 - 5 (3T J;+1 ) ln 2 ' 5 (3 In 3T jl+l ) Sl In 2 3 Sl 

1=1 

< e- x c 2 (Ni... N r ) 1/2 ln s+2 - 5r (2 r+1 A^ 1 ...N r ) ln 25r 3 ln(6iVi ■ ■ • N r ) (50) 



x 



with 



c 2 = s p + 3 5r 2.5 s s r2 ~ r/2 si • • • s r In 2 3s x • • • In 2 3s r 



It is easy to see that c 2 > max(p, 3000). Hence, if iVi • • • N r < max(p, 3000), 
then (!50j) is also true. Theorem 2 is proved. □ 
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